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humidity.  This  pap«r  will  fbcus  on  turbulenc*  proportios  and  turbulant 
heat  flux  measured  during  four  marine  fog  events.  The  fogs  were  charac- 
terized by  a slightly  unstable  boundary  layer  with  moderate  increases  of 
heat  during  the  fog.  An  ensemble  average  of  the  fog  events  showed 

that  (C^Wd  e reached  peak  values  isnediately  before  and  after  the  fog 
was  encountered.  An  analysis  of  the  dimensionless  temperature  structure 
^finietien'''^^^  normalized  by  the  square  of  the  temperature  gradient)  gave 
lower  values  of  tenperature  fluctuations  in  fog  than  in  clear  air. 
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Figure  Captions 


Figure  1 


Figure  2 


Figure  3 


Figure  4 


Figiire  5 


Figure  6 


Naval  Postgraduate  School  P/V  Acania  at  the  conclusion  of 
CEWCOM-76 . 


Nomalized  heat  conduction  coefficient  (eddy  dlffusivity) 
versus  Richardson  number. 

Variation  of  atmospheric  quantities  at  Z « 10  meters: 

a)  ten^rature  structure  function, 

b)  rate  of  dissipation  of  velocity  fluctuations,  e 

c)  wind  velocity,  U 

2 

Teiiq>erature  structure  function,  C,p  , at  Z > 10  meters 
versus  time  in  Monterey  Bay.  Solid  circles  are  unstable 
conditions  (12/76)  and  solid  triangles  are  near  neutral 
conditions  (7/76) . 
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Diurnal  variation  of  : 

a)  open  ocean 

b)  coastal  ocean 

R/V  Acwia  trades  off  Southern  California  during  fog  events 
for  CEWCOM-76; 


a) 

Sep  27-29 

(fog  #1) 

b) 

Oct  5-9 

(Fog  #'s  2,3,4) 

The  six  digit  number  indicates  the  day  of  the  month  (dd)  and 
time  of  day  (tttt)  in  the  code  - ddtttt. 

Sea  surface  temperature  contours  from  aircraft  XR  measurements. 
Courtesy  of  Backes  (1977)  and  Ralph  Mukson,  Airborne  Research 
Associates . 


Figure  7 


Figure 

8 

C ^ at 
T 

Z » 

10 

meters  versus  time  for  fog  events: 

a) 

Fog 

#1 

Sep  27-28 

b) 

Fog 

»2 

Oct  5 

c) 

Fog 

*3 

Oct  8 

d) 

Fog 

«4 

Oct  9 

Figure  9 

e at  Z 

» 10  meters  versus  time  for 

fog  events: 

a) 

Fog 

#1 

Sep  27-28 

b) 

Fog 

»2 

Oct  5 

c) 

Fog 

#3 

Oct  8 

d) 

Fog 

#4 

Oct  9 

Figiire 

10 

Sensible  heat 

flux,  W,  versus  time 

for  fog  events 

a) 

Fog 

#1 

and  Fog  #2 

b) 

Fog 

#3 

and  Fog  #4 

Figure  II  Radiosonde  measurements  of  teo^rature  (solid  line) 

and  relative  humidity  (dashed  line)  versus  altitude 
for  fog  events: 

a)  Fog  #1 

b)  Fog  #2 

c)  Fog  #3 

d)  Fog  #4 
Courtesy  of  NWSMET. 

Figure  12  Ensemble  average  of  turbulence  parameters  one  hour  before 

and  after  fog  is  encountered  and  one  hour  before  and  after 
fog  dissipates: 

a)  versus  time 


b) 


e/£  versus  time 


f 

i 

1 

2 

Figure  13  Height  dependence  o£  normalized  to  the  Z » 6.6  meter 

value  for  fog  (x's)  and  non-fog  (o's)  data.  Neutral 
equilibrium  yields  the  line  indicated,  non-neutral 

equilibrium  would  be  a less  steep  line  in  this  representa- 
tion. 

Fig\ire  14  Dimensionless  temperattire  structure  parameter  versus 

Richardson  number  (atmospheric  stability)  for  fog  (x's) 
and  non-fog  (o's)  data.  The  solid  line  Is  the  overland 
results  from  V^gaard,  et.  al.  (1971). 
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Atmospheric  Turbulence  Measurements  in  Meurine  Fog 
during  CEWCOM-76 

C.W.  Fairallf  K.L.  Davidson,  T.M.  Houlihan  and  G.E.  Schacher 

Environmental  Physics  Group 
Naval  Postgraduate  School 
Monterey,  California  93940 

Abstract 

Atmospheric  turbulence  parameters  e (turbulent  energy  dissipation 

rate) , and  (temperature  structure  function)  were  meaisured  in  the 

marine  boundary  layer  during  the  Cooperative  Experiment  for  West  Coast 

Oceanography  and  Meteorology  (CEWCOM-76) . The  atmospheric  stability 

was  characterized  by  measurements  of  profiles  of  mean  wind  velocity, 

temperature  and  relative  humidity.  This  paper  will  focus  on  turbulence 

properties  and  t\irbulent  heat  flux  measured  during  four  marine  fog 

events.  The  fogs  were  characterized  by  a slightly  unstable  boundairy 

layer  with  moderate  Increases  of  heat  flux  during  the  fog.  An  ensemble 

2 

average  of  the  fog  events  showed  that  and  e reached  peak  values 

immediately  before  and  after  the  fog  was  encountered.  An  analysis  of  the 

2 

dimensionless  tesperature  structure  function  (C^  normalized  by  the  square 
of  the  temperature  gradient)  gave  lower  values  of  temperature  fluctuations 


in  fog  tham  in  clear  aur. 


I.  Introduction 


During  September  and  October  of  1976  the  Naval  Postgraduate  School 
(NPS)  shipboard  turbulence  and  profile  systems  were  used  In  a marine  fog 
and  boundary  layer  research  cruise  off  Southern  California  aboard  the  NFS 
research  vessel  R/V  ACANIA  (Fig.  1).  Termed  the  Cooperative  Experiment  In 
West  Coast  Oceanography  and  Meteorology  (CEUC(M-76),  the  cruise  was  jointly 
organized  by  NPS  and  the  Naval  Electronics  Laboratory  Center  (NELC)  of  San 
Diego  In  collaboration  with  CALSPAN  Corporation  under  the  sponsorship  of 
NAVAIR.  Also  represented  on  board  were  the  Naval  Avionics  Facility, 
Indianapolis  (NAFI) , the  Naval  Research  Laboratory  (NRL) , the  California 
Air  Resources  Board  (CARS),  and  a Naval  Weather  Systems  Service  Mobile 
Environmental  Team  (NWSMET) . This  was  the  most  recent  of  a series  of  data 
gathering  cruises  forming  a major  part  of  a program  of  marine  fog  Investi- 
gation sponsored  by  Naval  Air  Systems  Command.  It  was  also  the  third  such 
cruise  to  Include  measurements  of  atmospheric  turbulence.  The  first  was 
conducted  in  1974  aboard  the  R/V  ACANIA  off  Northern  California  (Mack  et  al, 
1975);  the  second  was  conducted  In  1975  aboard  the  USNS  HATES  off  Novla 
Scotia  (Gathman  and  Larson,  1976). 

Marine  fog  Is  one  of  the  least  understood  and  most  complicated  atmos- 
pheric phenomena.  Mack  (1975)  has  Identified  at  least  five  different  fog 
formation  "scenarios"  over  the  ocean.  Although  considerable  fog  forecasting 
success  has  been  obtained  by  skilled  phenomenologlsts  employing  various 
fog  indices  (such  as  Inversion  height  and  air  sea  temperature  difference) 

^n  local  coastal  areas  (Leipper,  1948,  Backes,  1977),  accurate  marine  fog 
prediction  In  the  open  ocean  will  surely  benefit  from  applications  of  know- 
Leuge  from  the  synoptic  to  the  microscale.  On  the  microscale,  the  pro- 
duction of  fog  Is  governed  by  several  (Interacting)  factors:  aerosols. 


9 


t 


1 

radlaclon,  hydrostatic  stability,  temperature  inversion  characteristics, 
and  turbulent  transport  of  heat  and  water  vapor.  Turbulent  transport  was 
first  incorporated  into  fog  modeling  by  Rodhe  (1962).  More  modern  fog 
models  (Lee  and  Lin,  1975;  Barker,  1975)  employ  sophisticated  turbulent 
parameter  stability  and  scaling  laws  (see  section  II)  that  are  unverified 
in  Che  marine  boundary  layer  (Davidson  et  al,  1977).  The  main  thrust  of 
the  NFS  turbulence  group  during  CEWCOM-76  was  to  attack  these  problems  by: 

1)  providing  modelers  with  measurements  of  turbulence  and  mean  profiles 
in  £og,  2)  obtaining  data  for  evaluation  of  the  turbulence  scaling  laws 
in  the  marii  e boundary  layer,  3)  examining  some  of  the  characteristics  of 
turbulence  "background"  levels  for  non-fog  conditions  in  the  marine  boundary 
layer. 
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II.  Theory 

A.  Bouadary  layer  equations 

The  boundary  layer  is  chat  part  of  the  atmosphere  where  friction 
with  and  heating  by  the  surface  play  an  important  role  in  the  generation 
of  turbulence.  Near  Che  surface  Che  shear  stress  and  scalar  fluxes  are 
essentially  constant  with  height.  In  this  region  the  fluxes  can  be  rep- 
resented in  terms  of  scaling  parameters  (such  as  u^  and  T^)  that  are  in- 
dependent of  height.  We  shall  refer  to  this  layer  of  nearly  constant  shear 
stress  and  flux  as  the  surface  layer.  For  a complete  treatment  of  the 
surface  layer  equations  we  suggest  Lumley  and  Panofsky  (196A) , Buslnger 

(1971)  and  Kraus  (1972).  The  normalized  fluxes  of  momentum  (F  ),  heat 

m 


(F,  ) and  water  vapor  (F  ) are 
h q 


'n  ” ~ <'u'w'>  - u*^ 


\ - - <T’w’>  - u*  T* 


- (q'w'^  - u* 


(la) 

(lb) 

(lc) 


where  u'  and  w'  are  Che  fluctuating  horizontal  and  vertical  wind  velocities, 

respectively,  T'  is  the  fluctuating  potential  temperature  and  q'  is  the 

fluctuating  specific  humidity.  The  Reynolds  shear  stress,  T,  is  related  to 

2 

the  friction  velocity,  u^,  by  T =•  p u^  , where  P is  the  density  of  air.  In 
the  surface  layer  we  can  write  simplified  turbulent  energy  balance  equations 
using  production  terms  (proportional  to  the  vertical  gradient,  3x/3Z j and 
dissipation  terms  [e  ] as  follows: 


2 

velocity:  u^  3u/3z  - ^ u^  T^^  ■ e 
temperature:  u^  T^  3T/3z  ■ 
humidity:  u^^  q„  3q/3z  ■ e 


(2a) 

(2b) 

(2c) 
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Note  Chat  the  velocity  equation  (where  = e)  has  an  extra  pro- 
duction (or  suppression)  term  due  to  Che  bouyancy  of  air  (T^  is  the  vir- 
tual potential  temperature).  In  the  surface  layer,  the  appropriate  tur- 
bulent transport  scale  parameter  is  the  height  above  Che  surface,  z.  In 
this  regime,  turbulence  properties  can  be  scaled  through  the  similarity 
parameter  z/L,  where  L is  the  Monln-Obukhov  stability  length. 


T u 
<g  T * ’ 

V* 


(3) 


where  < ■ ,35  is  the  von  Karmon  constant.  We  can  now  define  the  mean  pro- 
file functions  <l>^  (z/L)  as  (Businger,et  al,  1971) 


(z/L) 


Kz  3x 
X*  3z 


(4) 


function  (Wyngaard  and  Cote,  1971) 

(5) 

Chen  Che  energy  equations  become 


(z/L)  - z/L  - (z/L) 

(6a) 

<j>T 

(z/L)  - E.J.  (z/L) 

(6b) 

(z/L)  - E (z/L) 

(6c) 

q q 

Present  estimates  of  and  are  given  in  Appendix  A. 

In  the  inertial  subrange,  turbulence  is  nearly  isotropic  and  we 
can  represent  the  Intensity  of  turbulence  fluctuations  as  a one-dimensional 
power  spectrum,  S^(k),  of  the  form  proposed  by  Kolmogorov 

S^(k)  - .25  k"^/^  (7) 


If  we  define  a dimensionless  dissipation 


e Ka 

E„  (z/L)  = 
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wlier«  k Is  the  wavenumber  and  Is  the  structure  function.  If  the 


quantity  x has  value  x(r)  at  position  r and  value  x(r  + d)  at  position 
r + d then 


C ^ x(r)  - x(r  + d)l^^  d 


•2/3 


(8) 


In  the  inertial  subrange  is  Independent  of  d.  Corrsln  (1951)  has  shown 


that  the  structure  function  can  be  related  to  the  dissipation  by 


C ^ - 4 6 e 

X X X 


(9) 


where  S is  a constant,  6-  ■ 6 * .81  and  6 ■ .52.  We  can  construct  a 

X 1 q u 


dimensionless  temperature  structure  parameter,  f(z/L),  by  combining  eqs. 
6 and  9 (Wyngaard,  et  al. , 1971) 


fl(z/L) 


(10) 


or,  using  eq.  4 


- OT/3z)2  2^/^  f2(2/L) 


(11) 


We  can  also  use  eq.  5 to  obtain  a similar  expression  for  £. 


e - — E (z/L) 

ica  “ 


(12) 


Since  the  Monln-Obukhov  length  is  difficult  to  measure,  we  normally  use 


an  alternate  stability  varlable-the  gradient  Richardson  number,  R^, 


(9V3z) 

*^1  “ T Ou/3z)2 


(13) 


J I 
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which  is  related  Co  z/L  by 


- (z/L)  ^^(z/L)/i<t>^(z/L))^  (14) 

2 

The  height  dependence  of  and  e can  be  obtained  from  eqs.  10  and  12;  for 

2 -2/3  -1 

near  neutral  conditions  (R^  - 0)  - z and  e - z . For  unstable 

2 -4/3 

conditions  (R^  « 0)  - z and  e «•  constant.  For  stable  conditions 

^ 0),  turbulent  eddies  are  increasingly  by  Increasing  stable 

scraciflcation  of  the  atmosphere.  Overland  measurements  Indicate  complete 
damping  of  turbulence  for  R^  > .21. 

B.  Flux  calculations 

There  are  four  methods  of  calculating  scalar  fluxes:  eddy  correlation, 
gradient,  turbulence,  and  bulk.  The  eddy  correlation  method  Is  a direct 
measurement  of  x'w*  In  equation  1.  This  requires  a stable  platform  so 
we  cannot  use  It  for  shipboard  measurements.  The  bulk  method  (Friehe,  1977) 
relates  the  flux  to  the  mean  wind  velocity  and  the  air-sea  temperature  or 
humidity  difference.  We  have  found  this  method  to  be  unsatisfactory  due 
to  difficulties  In  measuring  the  sea  surface  temperature.  The  gradient 
method  relates  the  vertical  gradient  of  the  scalar  to  the  flux  through  the 
eddy  dlffuslvlty,  For  instance,  the  normalized  heat  flux  Is 

Fh  - - K.J.  OT/3z)  (15) 

Unfortunately,  is  a function  of  stability, 

K.J,  • K z u*/4i^  (z/L)  (16) 

shown  In  Fig.  2.  Note  that  ILj,  - 0 for  Rj^  i *21  in  this  model.  The  tur- 

2 

bulence  method  Is  based  on  calculations  of  T^  from  and  e.  Combining 
eqs.  9 and  6b  we  find 


14 


(17) 


The  sensible  heat  flux  is  related  to  the  normalized  flux  by 

W » C F.  (19) 

h 

where  C - 1.3x10^  joule/ (‘C  m^)  is  the  heat  capacity  of  air.  Equations  15 
and  18  have  their  analogous  forms  for  calculation  of  water  vapor  flux.  Un- 
fortunately, Instrumental  limitations  and  problems  (see  Section  III)  have 

2 

prevented  us  from  measuring  C and  (9q/3z)  during  fog  events. 


I.  Instrumencaclon  and  Analysis 


Since  Che  NPS  work  represented  only  a part  of  Che  Cocal  shipboard 
effort  for  CEWCOM-76,  we  will  first  list  Che  key  measurements  made  by  Che 
ocher  groups: 

1.  CALSPAN  - aerosol  spectra  and  chemical  composition,  visibility 
and  profiles. 

2.  NRL  (R.  Jeck)  - Knollenberg  aerosol  spectrometer. 

3.  NAFI  (J.  Russell)  - microwave  refractometer. 

4.  CARB  - ozone,  carbon  monoxide  and  NO  . 

X 

5.  NWSMET  - twice  daily  radiosondes. 

The  NPS  Instrumentation  was  primarily  devoted  to  measurement  of  mean 

and  turbulence  profiles.  A comprehensive  description  of  this  system  Is 
the  subject  of  another  paper  (Houlihan,  et  al. , 1977).  The  profiles  were 
measured  at  four  fully  Instrumented  levels  plus  a sea  surface  temperature 
sensor.  The  level  heights  were  4.2,  6.6,  11.3  and  17.7  meters  above  the 
mean  waterline. 

Mean  temperature  was  measured  with  Hewlett-Packard  Model  RP2801A  quartz 
oscillator  thermometers,  accurate  to  .01  *C.  Mean  humidity  was  measured 
with  Hydrodynamic , Inc . Digital  II  type  15-1818  LiCl  Dunmore  sensors,  ac- 
curate to  about  3Z  RH.  These  humidity  sensors  have  very  slow  response 
above  about  95Z  RH  and  response  variations  from  sensor  to  sensor  make  them 
unreliable  under  fog  conditions.  For  this  reason,  we  do  not  have  fog 
event  humidity  profiles  of  sufficient  accuracy  to  calculate  the  humidity 
flux.  The  wind  velocity  was  measured  with  Thomthwaite  Model  101  cup 
anemometers,  accurate  to  about  3Z. 

Wind  speed  fluctuations  (for  e)  were  measured  with  Thermo-Systems,  Inc. 

Model  TSI  1054B  constant  temperature  anemometers  using  TSI  1210  probes  with 

L6 

_ J 


6 mil.  cylindrical  hoc-film  sensors.  The  hot  films  were  callbraced  on  a 

TSI  il2S  Calibrator  with  Che  cables  acCually  used  during  Che  shipboard 

2 

measuremencs.  Temperacure  fluccuatlons  (for  ) were  measured  with  GTE 
Sylvanla  ac  whmatatone  bridges  (3.0  kHz  carrier)  using  TSI  1210-P.8  micro 
thermal  probes  w|ch  2.5  y diameter  platinum  wire.  Ue  had  planned  Co  measure 
water  vapor  fluctuations  (for  ) ac  one  level  using  an  Electro-Magnetics 
Research  Corporation  Lyman-alpha  sensor,  however  this  instrument  failed 
early  in  the  cruise.  The  mean  data  was  continuously  averaged  and  period!- 
cally  logged  by  a microprocessor  designated  MIDAS  (Mlcroprogrammable 
Integrated  Data  Acquisition  System)  developed  at  NFS  and  adapted  for  the 
ACANIA  by  Plunkett  (1977).  The  fluctuating  data  was  amplified,  filtered 
and  recorded  on  a Honeywell  Model  5600  FM  tape  recorder. 

In  addition,  the  NFS  acoustic  radar  (Aeroenvlronment  Model  300)  was 
mounted  on  the  ship  and  provided  a contincous  monitor  of  the  Inversion 
height. 

The  analysis  of  the  turbulence  signals  was  quided  by  noise  consider- 
ations. The  combination  of  long  sensor  cables  (about  150  ft)  and  somewhat 
Impure  ship  generated  power  restricted  our  analysis  options — particularly 
for  the  hoc  film  systems.  The  velocity  signals  were  processed  with  Fourier 
power  spectrum  analyzers  to  yield  S^(f).  Using  Taylor's  hypothesis 
(k  " 2irf/u)  and  equations  7 and  9 we  find 

e - 2.67  ^ [f5/^  S^(f)]  (20) 

The  spectral  method  has  the  advantage  that  power  source  noise 

2 

(60,120,180  Hz,  etc.)  can  be  Ignored.  The  bulk  of  the  data  were  ob- 
tained from  paired  sensors  (separated  by  0.3  meters)  using  eq.  8.  In  some 

2 

cases  only  one  sensor  was  operable;  these  were  spectrum  analyzed  and 


was  calculated  using  eq.  7.  Being  based  upon  the  high  frequency  turbulence 
components,  both  of  these  techniques  are  relatively  unaffected  by  ship 


motions . 
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IV.  Results 

A.  General  Background 

Before  plunging  Into  the  fog  data,  It  is  of  Interest  to  examine 

the  variability  of  turbulence  parameters  In  the  marine  boundary  layer  to 

establish  a feeling  for  typical  "background"  behavior.  Nominal  variations 
2 

of  C.^  and  e during  a 24  hour  period  (Fig.  3a  and  Fig.  3b)  In  the  open  ocean 

are  relatively  smooth.  Since  c Is  proportional  to  u^,  we  expect  it  to  be 

1 

more  variable  than  . . A typical  dally  variation  of  wind  velocity  at 

2 

z > 10  meters  (u^q)  Is  shown  In  Fig.  3c.  Of  course,  Is  quite  sensitive 
to  changes  In  atmospheric  temperature  gradients  as  we  can  see  in  a com- 
parison of  measurements  In  Monterey  Bay  under  neutral  conditions  (7/76) 

2 

and  unstable  conditions  (12/76)  In  Fig.  4.  The  diurnal  variation  of 

overland  Is  dominated  by  solar  heating  of  the  surface  during  the  day  and 

2 

radiative  cooling  at  night.  This  leads  to  a maximum  of  In  the  afternoon 
and  minimum  around  sunrise  and  sunset.  Since  the  ocean  has  both  high  heat 
capacity  and  excellent  thermal  conductivity  (due  to  turbulent  mixing) , 
the  surface  ten^erature  Is  subject  to  only  a few  tenths  of  degrees  modu- 
lation by  solar  heating  cycles.  In  the  slightly  unstable  (R^  > - 0.1) 

boundary  layer  we  typically  encountered  In  the  open  ocean,  we  found  that 

2 

radiation  cooling  of  the  atmosphere  produced  a slight  maximum  of  during 
the  night  (Fig.  Sa).  In  coastal  areas,  influenced  by  land-sea  breeze  cycles, 
we  found  diurnal  variations  more  characteristic  (but  less  extreme)  of  those 
found  over  land  (Fig.  5b). 

B.  Fog  Data 

This  paper  will  concentrate  on  the  first  four  fog  events  encountered 
during  CEUCOM  76: 

Fog  n Sept.  27  1100  - 0800  (Sept.  28) 
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Fog  n 

Oct.  5 

0700  - 1100 

Fog  #3 

Oct.  8 

1000  - 1800 

Fog  i/4 

Oct.  9 

0800  - 1200 

Fog  event  data  are  provided  In  tables  In  Appendix  B.  The  location  of  the 

ACANIA  during  these  events  can  be  established  from  Fig.  6a  (Fog  #1)  and 

Fig.  6b  (Fog  #'s  2,3,  and  4).  Fog  #'s  1,3,  and  4 were  open  ocean  fogs, 

all  located  west  of  Point  Conception;  Fog  #2  was  a coastal  fog  located 

just  south  of  Palos  Verdes.  Sea  surface  temperature  contours  (Fig.  7)  for 

the  first  week  of  October  were  measured  by  Ralph  Markson  In  the  Airborne 

Research  Associate's  aircraft  using  a Barnes  PRT-5  Infrared  radiometer 

(Backes,  1977).  In  order  to  simplify  the  presentation  of  four  levels  of 

turbulence  data,  the  values  from  each  level  have  been  averaged  together 

after  translation  to  a 10  meter  equivalent  value  using  eqs.  10  and  12. 

The  10  meter  equivalent  averages  of  the  temperature  structure  function 
2 

(<C-  > ) for  the  four  fog  events  are  given  In  Fig.  8.  Note  the  conslder- 

10  2 

able  variability  of  compared  to  Fig.  3a  and  Fig.  4.  Figure  9 is  a 

similar  presentation  for  ^^^2.0*  some  cases  values  of  e are  not  avail- 

able during  fog  since  the  Impact  of  water  droplets  causes  a noise  spike 
by  cooling  the  film.  In  heavy  fog,  this  effect  totally  masks  the  real 

signal.  We  were  able  to  obtain  e values  during  Fog  #'s  1 and  4 by  using 

hot  wires  which,  due  to  their  much  smaller  size  (4.5  u diameter),  are 
less  affected.  The  heat  fluxes  for  the  fog  events  (Flg.s  10a  and  10b) 
were  calculated  from  the  10  meter  averages  using  eqs.  IS  and  18.  Rep- 
resentative radiosonde  dew  point  and  temperature  profiles  are  given  in 
Fig.  11. 

The  correlation  of  changes  in  turbulence  parameters  and  fog  events 
can  be  highlighted  by  preparing  "ensemble'^  averages  of  the  fog  encounters. 
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Each  £o(i  evenc  was  divided  Inco  four  one  hour  sections:  one  before  and 

one  after  the  fog  Is  first  encountered,  and  one  before  and  one  after  the 

fog  clears  (Fog  iH  was  excluded  from  this  analysis).  If  a given  fog  event 

was  less  than  two  hours  duration,  the  second  and  third  sections  were  shortened 

2 

to  1/2  the  fog  duration.  Note  that  both  and  e rise  to  a maximum  about 

30  minutes  before  fog  Is  encountered  (Figs.  12a  and  12b).  One  can  also 

2 ° 

examine  the  more  fundamental  behavior  of  In  terms  of  fog  and  non-fog 
comparisons  of  height  dependence  (eq.  10)  and  the  dimensionless  temperature 
structure  parameter  (DTSP),  fjCR^^),  from  eq.  11.  Although  the  height  de- 
pendence exhibits  no  obvious  dislmllarlty  (Fig.  13),  the  DTSP  curve  (Fig.  14) 
Implies  lover  values  of  for  a given  OT/3z)  during  fog  events.  This 
Is  consistent  with  the  lower  fluxes  obtained  during  fog  by  the  turbulence 
method  as  compared  to  the  gradient  method.  It  Is  quite  conceivable  that 
this  Is  an  Instrumental  effect,  caused  by  the  accumulation  of  liquid  water 
on  the  mlcrothermal  sensors.  Such  an  accumulation.  If  large  enough,  could 
reduce  the  frequency  response  of  the  sensors  sufficiently  to  attenuate  the 
measured  fluctuations.  Based  on  laboratory  measurements  of  frequency  re- 
sponse in  salt  coated  wires  (Falrall  & Schacher,  1977)  one  would  expect 

this  effect  to  be  small  until  the  water  film  was  on  the  order  of  30  microns 

2 

thick.  The  rapid  Increase  of  before  dissipation  of  fog  and  the  rapid 

decrease  after  dissipation  are  in  the  opposite  direction  of  an  expected 

water  film  effect.  Despite  the  scarcity  of  droplet  Impacts  observed  on 

the  130  micron  diameter  hot  films  (Implying  the  2.5  micron  diameter  micro- 

thermal  sensors  remained  very  dry)  during  Fog  #2,  a considerable  drop  In 
2 

was  measured. 

A 

Fog  event  data  from  other  investigators  can  be  found  in  Backes  (1977) 
and  Mack,  et  al.  (1977). 
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V.  Conclusions 


Despite  the  efforts  of  a number  of  scientists  during  a month  of  ship- 
board data  gathering,  unequivoceO.  statements  about  the  nature  of  marine 
fog  2u:e  still  difficult  to  come  by.  Although  it  might  be  wise  to  resist 
the  temptation  to  embark  into  arm  waving  arguments  about  the  reasons  for 
fog  formation,  it  is  interesting  to  note  that  the  three  open  ocean  fogs 
described  'were  found  in  the  vicinity  of  the  cold  water  intrusion  at  Point 
Conception  (upper  left  hand  corner  of  Pig.  6).  To  describe  these  fogs 
as  being  formed  by  "warm  air  blowing  over  cold  water"  is  probably  much 
too  optimistic,  since  all  were  characterized  by  unstable  surface  layers. 
Indeed,  the  vertical  heat  fluxes  often  dramatically  increased  during  the 
fog.  We  can  make  some  statements  about  the  nature  of  turbulence  in  fog, 
keeping  in  mind  that  they  are  subject  to  questions  about  the  influence 
of  water  droplets  on  sensors.  The  temperature  structure  function  (and 
therefore  the  turbulent  heat  transport)  was  subject  to  considerably  more 

variations  in  the  atmospheric  conditions  conducive  to  fog  formation  than 

2 

in  the  typical  non- fog  background.  An  ensemble  analysis  of  and  e 

showed  increased  turbulence  levels  just  before  and  just  after  the  fog  was 

encountered.  This  suggests  that  the  fog  is  often  super is^tosed  on  some 

large  scale  cellular  structure  in  the  boundary  layer  (on  the  order  of  10 

km  in  the  down  wind  direction) . When  analyzed  in  dimensionless  form 

2 

(normalized  by  temperature  gradient),  is  lower  in  fog  than  out  of  fog. 
It  is  possible  the  fog  droplets  are  retarding  the  ten^rature  fluctuations 
by  evaporation-condensation  response. 

In  general,  we  have  made  consideradsle  progress  towards  the  specific 
goals  enumerated  in  the  introduction.  Notably  in  the  provision  of  a 
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body  of  data  for  modelers  to  work  with.  However,  there  is  considerable 
room  for  improvement  (particularly  in  the  water  vapor  flux  capacity)  in 
future  work. 
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Appendix  A 

The  forms  of  cite  mean  profile  functions  are  taken  from  Buslnger,  et  al  (1971). 
<|>^(z/L)  » (1  - 15  z/L  < o 

(l>^^(z/L)  - (1  + 4.7  z/L)  z/L  > o 

4>t(z/L)  - .74  (1-9  z/L)'^^^  z/L  < o 

<J)^(z/L)  - (.74  + 4.7  z/L)  z/L  > o 

The  dimensionless  velocity  dissipation  function  is  from  Wyngaard  and  Cote 

(1571). 

E^(z/L)  - [1  + .5|z/l1^^^]^^^  z/L  < o 

E^(z/L)  - [1  + 2.5(z/L)^^^]^^^  z/L  > o 

The  dimensionless  temperature  structure  functions  is  from  Wyngaard,  et  al 
(1971). 

f3^(z/L)  - 4.9  (1-7  z/L)’^/^  z/L  < o 

fj^(z/L)  - 4.9  (1  + 2.75  z/L)  z/L  > o 
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